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Abstract 

In this note we investigate the effective superpotential of an A = 1 
SU{Nc) gauge theory with one adjoint chiral multiplet and Nf fun¬ 
damental chiral multiplets. We propose a matrix model prescription 
in which only matrix model diagrams with less than two boundaries 
contribute to the gauge theory effective superpotential. This prescrip¬ 
tion reproduces exactly the known gauge theory physics for all Nf 
and Nc- For Nf < Nc — 1 this is given by the Affleck-Dine-Seiberg 
superpotential. For Nf > Nc +I we present arguments leading to the 
conclusion that the dynamics of these theories is also reproduced by 
the matrix model. 


Introduction: In a recent series of papers iSi, Dijkgraaf and Vafa 
have proposed a perturbative method for computing the effective glueball 
superpotential of several classes of iV = 1 theories. This superpotential is 
essentially computed by summing over all the planar zero momentum Feyn¬ 
man diagrams of the theory. To better organize this computation in the case 
of helds in the adjoint representation of the gauge group, it is useful to ex¬ 
press it as a matrix path integral with the potential given by the tree level 
superpotential of the original theory. 

Although this duality was hrst obtained via a “string theory route” (build¬ 
ing on previous work in ^ , it is purely a held theoretic duality, and very 
recently it has been derived within a gauge theoretic framework for SU{Nc) 
theories with adjoints []^. Other related work on this duality has appeared 

in i 1,0,0,0, lii, g. 

One of the most natural extensions of this duality is to theories with 
helds in the fundamental representation of the gauge group. The theory 
under consideration is supersymmetric QCD with gauge group SU{Nc) and 
Nf havors coupled with a single chiral held in the adjoint representation. 

One way to proceed is to represent the Nf helds in the fundamental 
representation as a Me x Mf matrix while keeping Mf/Mc = Nf/Nc, and 
extend the matrix integral to contain such objects. Perturbing the resulting 
potential by mass terms for these helds allows one to integrate them out hrst 
and obtain an ehective superpotential for the remaining adjoint helds. This 
method has been proposed in [g, and has yet to be tested against gauge 
theory predictions. 

Another method [0] is based on the same integral, but uses the well- 
known fact that a quark loop in a large Nc graph corresponds to adding 
a boundary to the corresponding Riemann surface. If one has Nf havors, 
graphs with multiple boundaries are weighted by powers of It is there¬ 
fore possible to argue that, in the limit of small Nf/Nc, the matrix integral 
is dominated by two terms, one without quark loops, which is the typical 
Veneziano-Yankielowicz term and is proportional to the rank Nc of the gauge 
group, and the other coming from the diagrams with a single boundary. In 
this limit the contributions of multiple boundaries are suppressed. However, 
for Nf/Nc of order unity it is not a priori clear in this prescription that the 
diagrams with a higher number of boundaries are negligible. 

This method was successfully used to compute in the matrix theory the 
effective superpotential of an SU{2) gauge theory with one flavor |g. Very 
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recently this analysis was extended to include theories with gauge group 
SU{Nc) and one flavor |^. The superpotential evaluated at one of its ex¬ 


trema was presented as a series in the bare Yukawa coupling and was shown 
to agree up to 7’th order in the expansion parameter with the one obtained in 
the gauge theory after integrating out the adjoint held, adding the Affleck- 
Dine-Seiberg (ADS) nonperturbative contribution and integrating out the 
quark fields. 

The puzzle that triggered our analysis was the discrepancy between the 
apparent possibility of considering a larger number of boundaries in the ma¬ 
trix model analysis and the impressive agreement between gauge theory and 
the 1-boundary superpotentiah Indeed, in the analysis of [IT diagrams with a 
higher number of boundaries are suppressed only by powers of Nf/Nc = 1/2. 
Since the gauge theory result is exact, and the agreement is exact up to fac¬ 
tors of the order of 1/2*^ where n is a large number, one cannot but suspect 
that this too good an agreement is due to the fact that all matrix model 
diagrams with more than a single boundary do not contribute at all. 

We are therefore led to formulate a prescription in which only diagrams 
with less than two boundaries are relevant. This is accomplished by repre¬ 
senting the gauge adjoint fields as M x M matrices, and the Nf fundamental 
fields as M X Nf natrices, and taking the large M limit. In the case when no 
fields in the fundamental representation are present, this reduces to the orig¬ 
inal Dijkgraaf-Vafa proposal. Our proposal is very similar in spirit to the one 
TT[] ; the crucial difference is the automatic suppression of contributions of 


m 

matrix diagrams with more than two boundaries. 

The main puzzle raised by our prescription is that the matrix model does 
not distinguish between Nf < Nc and Nf > Nc while the gauge theory 
physics changes drastically. 

In this paper we first show that the effective superpotential obtained in 
the gauge theory by integrating out the adjoint held, adding the ADS super¬ 
potential, further perturbing by quark mass terms, and then integrating out 
the quarks is in exact functional agreement with the effective superpotential 
given by the matrix model graphs with less than two boundaries, after inte¬ 
grating out the glueball superfleld S. Thus, the gauge theory superpotential 
(with the ADS component included) and the superpotential obtained in the 
matrix model are obtained from the same effective superpotential W ^(A) 
by integrating in different fields. 

We then argue that, for Nf > Nc, the gauge theory superpotential eval- 
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uated at its extrema is also reproduced by the matrix model. 

The motivation behind our proposal was the impressive agreement be¬ 
tween gauge theory and the matrix superpotential due to diagrams with less 
than two boundaries. However, our computations do not exclude the (un¬ 
likely) possibility that matrix diagrams with more than two boundaries are 
absent only in the particular case we consider. We will comment on alterna¬ 
tives at the end of our paper. 


Gauge theory: It was shown by Seiberg a long time ago that the re¬ 
quirement of holomorphy as well as the unbroken symmetries can be powerful 
allies for hnding exact results regarding effective superpotentials in = 1 
supersymmetric theories. The simplest example of this sort is supersymmet¬ 
ric QCD. This theory contains an SU{Nc) vector multiplet as well as Nf 
pairs of quark helds, Q, Q, with Q in the antifundamental of SU{Nc) and 
the fundamental of the SU{Nf) flavor group, and Q in the fundamental of 
SU{Nc) and the antifundamental of a different SU{Nf) flavor group. 

The theory has no tree-level superpotential; the Lagrangian is thus: 

C = J d^eTr[Qe^Q + Qe-^Q] + J d^OTrlW^Wa] (1) 


where the hrst and second traces are in flavor space while the last is in color 
space. 

There are several different ways of assigning charges to the various helds 
and coupling constants present in this theory. One of them, which follows 
by requiring that the anomaly cancellation condition is satished by the phys¬ 
ical helds alone, yields the following representations for the various helds, 
coupling constants and dynamical generated scale: 



SU{Nf) 

X SU(N,) 

X f/(l), 

X I'll). 

X U{1)r 

Q: 

Nf 

1 

1 

0 

Nf-N^ 

Nf 

Nf-Nc 

Nf 

Q: 

1 

Nf 

0 

1 

> 

CO 

1 

1 

1 

1 

1 

0 


Standard nonrenormalization theorems imply that there is no superpo¬ 
tential generated perturbatively. Nonperturbatively however, nonrenormal¬ 
ization theorems fail and a superpotential is generated ||I8|, IS. Its form, up 


to numerical coefficients, is hxed completely by symmetries. For Nf < N^. 
this can be argued in 3 steps. 
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1) First, one notices that the moduli space can be described by the gauge 
invariant meson helds Xd = Qi°'Qj- Requiring invariance under both flavor 
symmetry groups implies that the effective superpotential can be a function 
only of det X. 

2) The requirement that the R-charge of the superpotential is 2 implies 

that the exponent of detX is ■ 3 N^-Nf 

3) Finally, on dimensional grounds, one has to add A ^ where A is 
the dynamically generated scale of the theory. 

The fact that the dynamical scale of the theory does not appear under a 
logarithm implies indeed that this superpotential is nonperturbatively gen¬ 
erated. The coefficient of the resulting term can be explicitly computed for 
Nc — Nf = 1 ([jT9[ for the SU{2) gauge group and for SU{Nc)) and is 
found to be equal to unity. For — Nf > 1 the coefficient can be obtained 
by adding mass terms to the appropriate number of flavors and integrating 
them out. The result is the usual ADS superpotential [18|: 


Wdyn = (X, - Nf) 


' Jy3Nc—Nf \ Nc-Nj: 


detX 


(3) 


A slight deformation of the theory we are discussing is the inclusion of a 
tree level superpotential: 


^aTr[QQQQ] = ^aTr[X^] (4) 

For consistency with the various symmetries of the original theory, the charge 
assignments for this new coupling constant are: 

SU(N,) X SuTn,) X 17(1), X (7(1). X (/(1)h 

„ . 1 1 o o 


This superpotential can be thought of as arising from integrating out a mas¬ 
sive adjoint held with a cubic coupling with the quarks. 

One can show that in this theory the only superpotential that can be 
generated is again the ADS one. Furthermore, adding mass terms to some of 
the quarks and integrating them out produces a new effective action which 
has the same form as the original one. This is achieved by keeping the 
coupling constant a fixed while taking the quark mass to infinity and keeping 
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a certain combination of this mass and dynamically generated scale fixed. 
This procedure also yields the change in the dynamical scale to be: 




( 6 ) 


where m is the mass of the quarks which were integrated out. 


The central object W 


crit 


that contains all the information about the 


dynamics of the theory can be obtained by integrating out all the helds. 
The resulting function depends on the scale Aq Q as well as on the coupling 
constant a and the masses of the quark £elds[[ Knowledge of this function 
allows one to reconstruct the full superpotential via Legendre transforms |TB |. 

Without loss of generality we can integrate out all the quark fields at 
the same time, by introducing a mass matrix m proportional to the identity 
matrix. It is certainly possible to introduce hierarchical masses (and/or non¬ 
diagonal mass matrices) and integrate out one quark field at a time. However, 
the hnal result will be expressed only in terms of the scale Aq; the expression 
in terms of more general mass matrices can be trivially restored. 

It is customary to write the superpotential in terms of the meson field 
X = QQ. Then, the superpotential deformed by mass terms is: 


W^ff = mTr[X] - -Tr[X^] + (W - Nf) 


A 


3N^-Nf 

_ 

detX 


(7) 


where is the scale of the theory with Nf light fields. Since we assumed 
that the mass matrix is proportional to the identity matrix, it follows that, 
at the minimum of HA//, X will have the same property. Thus, we can write 
from the outset 

X = xIn^ . (8) 

By simple inspection of the superpotential above it is clear that explicitly 
finding its extrema for general numbers of colors and flavors is not an easy 
operation. The approach we take is to use the equation of motion for x to 
cast the superpotential into a simple form in which one still has to replace x 
by the solution of some algebraic equation. 

’^The index denotes the fact that there are no more light fields in the theory. 

^If one interprets the quartic tree level superpotential as arising from integrating out 

an adjoint field, then the coupling constant a is a function of the mass of the adjoint field 

~ 2 

and the strength g of the trilinear coupling Tr[Q(l)Q] ; a = ^. 
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The equation for x is simply: 


3Nc-Nf 


A 


Nf 


m 


ax 


= 0 


XNC-Nf 

It is convenient to introduce new dimensionless variables 


(9) 


y = 


m 




X , p = 


m 


m 


^ 

m?M ° ’ 


k = 


Nr 


Nc-Nf 


, ( 10 ) 


where Aq is the dynamical scale of the theory with no light fields^. In terms 
of these variables the equation of motion for x becomes: 


Pyk + l_yk^^^Q , ( 12 ) 

or equivalently 

y-^ = l-Py , (13) 

since ?/ = 0 is not a solution of equation (^). 

Now we use these equations to eliminate the quadratic term in Wejj- A 
small amount of algebra leads to the superpotential evaluated at its minimum 
in terms of the newly introduced variables y and k: 


W 


1 

[ k + 1 1 ] 

= -NfAr 
crit 2 “ 

^ k — 1 


(14) 


where again y is a solution of equation ([T^) . 

This is the form that will be compared with the matrix model predictions. 
In the next section, using our modified extension of the Dijkgraaf-Vafa pre¬ 
scription to include fields in the fundamental representation, we will compute 
the value of the superpotential at its critical points and recover equations (0 
and (^) for all Nc and Nf. We will then present a gauge theoretic argument 
that (B and ([T3| ) are also true for Nf > Nc. Q 

^ Indeed, by taking iVj- = 0 in equation (^) and using the definition of k it is easy to 
see that 


Aq = m 


\2k+l 

m 


( 11 ) 


"'^The case = Nf is problematic since the ADS superpotential does not admit a 
continuation to this point. 
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The matrix model: The prescription of Dijkgraaf and Vafa instructs 
that to compute the effective superpotential for the gaugino condensate, we 
have to compute the planar partition function for the matrix model with a 
potential which is the tree level superpotential of the N = 1 theory we are 
interested in. The original arguments covered theories which had helds in 
bi-fundamental representations of the gauge group. 

Applied to a superpotential with a single critical point, this proposal 
yields the effective superpotential for the glueball superfield S by the follow¬ 
ing three steps: 

First, one computes the contribution to the free energy due to planar dia¬ 
grams, JFq. This is accomplished by formally replacing the gauge theory fields 
with M X M matrices in the matrix model potential. Then one computes 
the path integral in the limit M —^ cx). 

The second step is to identify the’t Hooft coupling in the matrix model 
with the gauge theory glueball superheld. At this stage JFq becomes a function 
of S only. 

The third and last step if to construct the gauge theory effective super¬ 
potential as: 

BT 

Wdv = N,^ + tS . (15) 


The only gauge theory ingredients entering this relation are - the number 
of colors, and r - the bare coupling constant. There is no relation between Nc 
and the dimension of the matrices used in the matrix model computations; 
the hrst one is a parameter while the second one was identihed with the 
glueball superheld. 

In the case of a trivial superpotential the only contribution to JFg comes 
from the normalization of the matrix path integral by the volume of the 
gauge group. In this case Wdv reproduces the Veneziano-Yankielowicz su¬ 
perpotential. 

Attempts to extend this prescription to include helds in the fundamental 
representation of the gauge group were formulated in |^, lOT. The idea is 


again to use the tree level superpotential as potential for the matrix model. In 
one proposal both the number of havors Nf and the number of colors Nc are 
“promoted” to matrix model variables Mf and Me such that Nj/Nc = Mf/Me 
[TSl . In the other proposal none of them is |^. In both cases, in the limit 


Nf Nc one cannot perform an expansion in the number of boundaries in 
the matrix model, which makes comparison with gauge theory difficult. 
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Here we propose that the matrix model variables are M x M matrices if 
the corresponding gauge theory helds are in the adjoint representation and 
M X Nf matrices if there are Nf helds in the fundamental representation 
in the gauge theory. Then, as in the original DV prescription, we identify 
the gauge theory glueball superheld with the’t Hooft coupling of the matrix 
model, i.e. S = g‘^M. As put forward in |]T3, the extended DV superpotential 


is the sum of the original one and extra contributions coming from matrix 
diagrams with boundaries: 


Wnv = K 


dJ^o 

dS 


tS J~i 


boundaries 


( 16 ) 


If one recalls that S was identihed with the dimension of the matrices 
and examines the M dependence of the various terms in the equation (|T^), it 
becomes clear that only diagrams with one boundary contribute in the planar 


(large M) limit. Indeed, equation (IT) can be organized as an expansion in 
[Nf/M). Then, due to the derivative with respect to S, the hrst term and 
the contribution of diagrams with one boundary will be of the same order in 
an 1/M expansion. All diagrams with two or more boundaries are suppressed 
by factors of 1/M and do not contribute in the large M limit. 

As in the original DV proposal, the part of jFg arising from the volume 
of the gauge group yields the Veneziano-Yankielowicz superpotential. This 
part is not modihed by the inclusion of helds transforming in representations 
of havor symmetry groups because these symmetries are only global. 

For supersymmetric QCD with an adjoint held 0 this leads to II. 


16 


= 


D.#>D(3D(3e“ 


(17) 


with a = ^, as in the gauge theory discussion. The path integral above can 


be computed in both of its forms, using the analysis of |^. Using the hrst 
form above and the analysis in we can show that the contribution to the 
free energy of the matrix model in the large M limit is: 


•^x=i = -NfS 


1 

2^ 4aS' 


In 


- + -yr 
2 2 


AaS 


( 18 ) 


where a = -%. Then, according to our extended DV prescription, the super¬ 
potential is: 


Wdv = N,S 


(19) 





























N,S{l-\n-^^)-N,S 


i -4aS - 1) -lni(l + \/l - 4aS) 


for all Nf and Nc. 

As in the gange theory case, we proceed by integrating ont the massive 
fields. Various terms in the equation above imply that the gaugino helds 
acquired a mass and thus can be integrated out. Despite the complicated 
form of the superpotential, its extrema are given by a very simple equation: 

N,\n-^ = Nf\n^{l + Vl - AaS) . (20) 

The hrst step is to use this equation to eliminate the logarithms in Weff- 
Then, the critical point equation can be cast in a more useful form: 


which in turn leads to the following expression for the superpotential at the 
critical points: 


W 


crit 


2a 


k + 1 
k-1 


(aS) 



+ 1 


( 22 ) 


This equation does not appear similar to the corresponding one on the 
gauge theory side. It is nevertheless possible to relate them more closely. To 
this end we massage equation A small amount of algebra combined 

with the observation that S' = 0 is not a solution of equation (0 , casts it in 
the following form: 


1 - 




(23) 


Introducing the notation 


w = 


{asy 


equation (|23|) becomes: 


w = 1 — (aA'^)t(; 


(24) 

(25) 
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which is the same as equation ([TB|). 

In terms of this new variable the superpotential evaluated at its minimum 
is given by: 


W 


crit 



k+l 1 

I—+ '«^ 
k — ^ 


(26) 


where w must be replaced by a solution of equation (p5|). 

Since w in the matrix model result as well as y on the gauge theory side 
are dummy variables, the results of the two computations agree provided 
that the matrix model scale A is identihed with the gauge theory scale in the 
conhning vacua as: 


= Aq = (detm)^A^ , (27) 

which is consistent with both theories not having massless helds. Further¬ 
more, since the analysis above goes through ii N^ this seems to suggest 

that the same should hold in the gauge theory as well. 


Nf > Nci It is clear from equation (]I^, and it was also pointed out in 
|1T6| that the extended DV superpotential is insensitive to the range of Nc 
and Nf. It was suggested in [T^ that this “problem” might be solved by 
including contributions from additional boundaries, which in the Nf ^ Nc 
limit would become important. The proposal we put forward in this paper 
implies that such contributions are inexistent. 

We have shown above that the 1-boundary matrix integral reproduces 
the expected gauge theory answer for all Nf and Nc for which the ADS 
superpotential holds. Moreover, this range can be extended to include also 
the case Nf = Nc + 1. Indeed, without introducing baryonic sources, the 
superpotential in this case is just the continuation of the ADS one 
this number of flavors. 


22 to 


The IT 


crit 


obtained in the matrix model always looks as if one had blindly 


started in gauge theory from an ADS-like superpotential, even for Nf > Nc- 
Our proposal implies that, even though for A^j > A(, -|- 1 no superpotential 
is generated in the beginning, once all quark helds are integrated out the 
superpotential evaluated at the critical point is still given by equations (p6[) 
and (^). At hrst glance this might seem problematic. It is nevertheless 
possible to argue that there is no contradiction with the held theory analysis 

0 . 

The argument is a slight generalization of the analysis in 


23 and is 


based on Seiberg’s duality together with the possibility of freely passing from 
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electric variables to magnetic variables in the path integral. 

In gauge theory, one can deform the theory by adding mass terms for all 
quark helds and then integrate them out. For Nj > Nc +2 the gauge theory 
is strongly coupled and the correct description is given by its Seiberg dual. 
In this description we could add masses to the magnetic quarks and integrate 
them out, until we reach the theory with a completely broken gauge group. 
The superpotential for this theory is written in terms of the helds which are 
dual to the electric meson helds and, in the absence of baryonic sources, is 
just the ADS superpotential. At this point we dualize back to the electric 
theory which is now weakly coupled. 

The results of the gauge theory analysis for Nf = Nc + I, given by (0 
and (0), can now be applied without reservations. Since the matching of 
scales when helds are integrated out is hxed by the renormalization group 
equations, it follows that the ehective scale when all quarks are integrated 
out is given by equation with the initial Nf and Nc [0, p3[| . 

Thus, by integrating out all helds in the fundamental representation. 


we obtain the same W 


crit 


as if we blindly started with an electric ADS- 


like superpotential for Nf > Nc + 2 and integrated out all quark helds. 
This explains the apparently unnatural agreement between the matrix model 
results and the introduction of an ADS-like superpotential for Nf > Nc + 2. 

Discussion: In this paper we formulated an extension of the Dijkgraaf- 
Vafa proposal which includes helds in the fundamental representation. Basi¬ 
cally, this proposal states that only matrix model diagrams with less that two 
boundaries contribute to the “gauge theory - matrix model” duality. We have 
tested this proposal by analyzing in detail supersymmetric QCD coupled 
to an adjoint chiral multiplet, both from the gauge theory and the matrix 
model perspective. We have found that, in the range of parameters where 
the Afheck-Dine-Seiberg superpotential is valid, it is reproduced exactly by 
the matrix model computation. Using Seiberg’s duality we have then argued 
that this agreement persists even for Nf > Nc, when the gauge theory is 
strongly coupled. Therefore the gauge theory and the matrix model analysis 
described earlier in this note agree for all values of Nc and Nf. 

Since for theories with Nf > Nc the baryons acquire a nonvanishing 
expectation value, it would be interesting to repeat the analysis by including 
baryonic sources both in the matrix model and in the gauge theory. This 
would extend the number of arguments of W and would provide a further 

crit 
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test of the extension of the DV prescription proposed in this note. 

It wonld also be interesting to see how onr proposal wonld be implemented 
when the recent held-theoretic proof of the initial DV proposal is extended 
to inclnde helds in the fnndamental representation. Another worthwhile en¬ 
deavor wonld be hnding the gange theory interpretation of matrix model 
diagrams with more than one bonndary. According to onr proposal they are 
snppressed like the nonplanar diagrams in the original matrix model; it seems 
however likely that they do not correspond to gravitational corrections, given 
the difficnlties with the geometric engineering of snch theories. 

As we pointed ont in the Introdnction, onr compntations do not exclnde 
the (nnlikely) possibility that matrix diagrams with more than two bonnd- 
aries are only absent in the particnlar case we consider. We can imagine 
several possible reasons for this: 

1) The extension of the DV proposal to inclnde fnndamental helds is still 
snch that only diagrams with a single bonndary contribnte. Nevertheless 
this conld be dne to mnltiple bonndaries being considered as disconnected 
diagrams, and thns not contribnting to the free energy. 

2) Diagrams with more than one bonndary simply vanish. If this were 
the case, this would give a highly nontrivial prediction for the matrix model. 

Whatever the cause, we believe that the impressive agreement between 
the gauge theory and matrix model results strongly supports the fact that 
matrix model diagrams with more than one boundary do not contribute to 
the duality, and we look forward to seeing how this phenomenon will emerge 
within a purely gauge theoretic framework. 
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